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TESTING THE MASSIVE DISK SCENARIO FOR IRAS 18089-1732 



in 
o 
o 

(N 
u 

< 

(N 



> 
oo 

o 
in 
o 

Or 

6 



13 



H. Beuther, Q. Zhang, T.K. Sridharan, Y. Chen 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 

hbeuther@cfa.harvard.edu 
Accepted by the Astrophysical Journal 

ABSTRACT 

Investigating in more detail the previously suggested massive disk scenario for the High-Mass Pro- 
tostellar Object IRAS 18089-1732, we observed the source in the 860 /mi band with the Submillimeter 
Array in various spectral lines and the submm continuum emission at (sub-)arcsecond spatial resolution. 
Fifty spectral lines from eighteen different species spanning upper level energy states between 17 and 
747 K were detected. One of the assumed best tracers for massive disks, CH 3 CN, is optically thick and 
does not allow a further disk investigation. However, the complex molecule HCOOCH3 appears optically 
thin and exhibits a velocity shift across the central core perpendicular to the emanating outflow. This 
signature is comparable to well-known low-mass disks and confirms the detection of a massive rotating 
structure likely associated with the central accretion disk. Assuming equilibrium between centrifugal 
and gravitational force, the estimated mass for this rotating structure is 16/(sin 2 (i)) M Q (with i the 
unknown inclination angle), of the same order as the gas mass derived from the continuum emission. 
Therefore, in contrast to low-mass disks, a considerable amount of the central gas mass is associated 
with the rotation, implying that massive disks may not be Keplerian. A temperature estimate based 
on the CH3CN(19-18) K-ladder results in ~ 350 K, thus a hot core has already formed in this region. 
An analysis of the submm continuum emission reveals even at this high spatial resolution only one iso- 
lated massive dust core without any detectable companions down to mass limits between 0.2 and 3 M Q 
(depending on the assumed temperature). Potential implications for the massive cluster formation are 
discussed. The radial intensity distribution of the core is highly non-symmetric, outlining the difficulties 
of density structure investigations based on radial averaging. 

Subject headings: star: formation - submillimeter - techniques: interferometric - accretion, accretion 
disks - astrochemistry - ISM: individual(IRAS 18089-1732) 



1. INTRODUCTION 

The phenomenon of massive protostellar disks is one of 
the central questions in high-mass star formation. Indirect 
indicators for the existence of massive disks has been ac- 
cumulated over recent years but unambiguous evidence for 
massive disks is still missing. High-spatial-resolution ob- 
servations of massive molecular outflows revealed that they 
could be as collimated as their low-mass counterparts (e.g., 
Beuther et al. 2002a; Gibb et al. 2003; Garay et al. 2004; 
Davis et al. 2004) indicating some underlying protostar- 
disk interaction as the driving mechanism. Furthermore, 
there have been various approaches to study massive disks 
in a more direct way. Infrared and mm continuum ob- 
servations find interesting disk-candidates (e.g., Shepherd 
et al. 2000; Chini et al. 2004), but they lack the necessary 
velocity information to provide kinematic proof of poten- 
tial accretion disks. Contrary to this, Class II CH3OH and 
H2O maser emission found disk- like kinematic signatures 
in a few cases (e.g., Torrefies et al. 1996; Pestalozzi et al. 
2004), but since maser emission is very selective and re- 
quires special conditions such studies do not allow a more 
statistical investigation of massive disk properties. The 
most promising approach to find and investigate disks in 
high-mass star formation appears to be the observation 
of thermal line emission which should be sensitive to the 
gas properties and kinematics of the disks. Various in- 
vestigations have been directed this way: one of the first 
massive disk studies were the CH3CN observations toward 



IRAS 20126+4104 which found a velocity gradient perpen- 
dicular to the outflow (Cesaroni et al. 1997, 1999). More 
recently, Beltran et al. (2004) also used CH 3 CN observa- 
tions to analyze potential rotating structures toward two 
even more massive high-mass star-forming regions, and 
again they found velocity gradients perpendicular to ob- 
served molecular outflows. A different approach was fol- 
lowed by Zhang et al. (1998, 2002) who investigated the 
disks of IRAS 20126+4104 and AFGL5142 by the means of 
NH3 observations. While these approaches are promising, 
there is still no unambiguous evidence for disks in high- 
mass star formation. Even more troubling is the fact that 
we still do not know which is the best molecular transition 
to study massive disks, CH3CN seems to be a promising 
candidate and NH3 works in some cases as well, but we 
also know a number of sources where the NH3 observations 
apparently do not trace the disks (Zhang, priv. comm.). 
One of the challenges of the coming years is to establish 
a reliable procedure to find massive disks and study their 
properties, which is important to constrain any theory of 
high-mass star formation. 

With this problem in mind, we began high-spatial- 
resolution observations of the HMPO IRAS 18089-1732 
in various molecular lines with the Submillimeter Array 
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(SMA 1 ) last year (Beuther et al. 2004b, a). This source at 
a distance of ~ 3.6 kpc with a luminosity of ~ 10 4 5 L© is 
part of a sample of 69 HMPOs intensively studied over 
recent years by Sridharan et al. (2002); Beuther et al. 
(2002b, c,d). The single-dish millimeter continuum obser- 
vations reveal a gas mass of ~ 2000 M Q , and H 2 and 
Class II CH 3 OH maser emission as well as a weak cm 
continuum source (0.9 mJy) has been detected at the core 
center. Previous SMA observations found an outflow in 
north-south direction and HCOOCH3 emission centered 
at the mm continuum peak exhibiting a velocity gradi- 
ent approximately perpendicular to the outflow (Beuther 
et al. 2004a). The latter observations were interpreted as 
indicative of a possible rotating disk-like structure around 
the central protostar which could be influenced by the out- 
flow as well as infall. The SiO outflow emission is slightly 
offset from the mm continuum peak (~ 1"), allowing to 
speculate about a potential secondary source undetected 
so far. Observations of a series of other molecular transi- 
tions revealed that many of them trace to some degree the 
outflow as well, and only a few remain confined to the core 
center, CH3OH being one of the latter species (Beuther 
et al. 2004b). However, some of the typical disk-tracers 
like CH3CN were not covered by the spectral setup. Fur- 
thermore, the spatial resolution was marginal for our pur- 
poses (2.7" x 1.7" corresponding to approximately 7000 AU 
at the given distance). Therefore, we followed up on the 
previous observations via higher spatial resolution (< I") 
observation covering a CH3CN series as well. These obser- 
vations are intended to get a better understanding of the 
potential massive accretion disk as well as to investigate 
which molecular lines are best suited for future massive 
disk studies. 

2. OBSERVATIONS 

We observed the HMPO IRAS 18089-1732 in 3 nights 
during the summer of 2004 with the SMA in two configu- 
rations covering unprojectcd baseline lengths between 20 
and 340 m. The phase center was R.A. 18:11:51. 4/Dec. - 
17:31:28.5 (J2000) and the tuning frequency 349.45 GHz 
in the upper sideband («i sr = —33.8). Two out of the 
three nights the weather was good with a zenith opac- 
ity of 0.18 at 348 GHz - scaled from the 225 GHz opacity 
measurement of the NRAO tipping radiometer operated 
by the CSO . In the third night, the weather was worse 
with t(348 GHz) ~ 0.38, the latter data were weighted 
considerably lower in the data reduction. The spectral res- 
olution over the whole bandpass of 4 GHz (upper + lower 
sideband) was 0.8125 kHz corresponding to a velocity res- 
olution of ~ 0.7km/s. The initial flagging and calibration 
was done in the IDL superset MIR originally developed 
for the Owens Valley Radio Observatory (Scoville et al. 
1993) and adapted for the SMA 2 . The imaging and data 
analysis was conducted in MIRIAD. The passband cali- 
bration was derived from a variety of different sources for 
each track including Jupiter, Neptune, Uranus, Callisto, 
Mars and 3c279 (for the largest baselines). Flux calibra- 
tion was conducted with observations of Uranus, Callisto 
and 3c279 (again for the largest baselines). The flux den- 

1 The Submillimctcr Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astron- 
omy and Astrophysics, and is funded by the Smithsonian Institution and the Academia Sinica. 

2 The MIR cookbook by Chunhua Qi can be found at http://cfa-www.harvard.edu/~cqi/mircook.html. 



sity scale is estimated to be accurate within 15%. For the 
phase and amplitude calibration we used regularly inter- 
leaved observations of the quasar NRAO530 (about 10° 
from the target source) with an approximate flux of 1.5 Jy 
at the given frequency. We applied different weightings 
for the continuum and line data resulting in synthesized 
beams for the continuum of 1.1" x 0.5" (P. A. 1.3°) and for 
the lines of 1.3" x 0.9" (P.A. 2.2°). The la rms of the 
continuum image is 35 mJy. 

3. RESULTS 

3.1. 860\im continuum emission 

Figure 1 shows the continuum image with 1.1" x 0.5". 
In spite of some structure, it remains a single source and 
does not split up into multiples. Two corresponding am- 
plitude versus uv-distance plots are shown in Figure 2. 
While the shortest baselines show a halo component, the 
range between 50 and 90 kA appears rather flat consistent 
with the previously modeled point source (Su et al. 2004). 
However, from 90 kA onward the flux densities decrease 
further, indicating additional sub-structure on the corre- 
sponding spatial scales. Beyond 200 kA, the S/N ratio is 
low, and in the 5kA bin plot (Fig. 2 left) one cannot dif- 
ferentiate the emission from statistical noise. Averaging 
the data into larger 40 kA bins to increase the S/N ratio 
(Fig. 2 right), we can discern emission above the noise. 
The large baseline data are consistent with a point source 
of approximate flux density ^0.3 Jy. The integrated flux 
density of the shorter baseline halo component can be fit- 
ted with a Gaussian to ~ 6.0 Jy. 

As a continuation of our radial profile study toward 
this source with single-dish 1.2 mm continuum observa- 
tions (with the MAMBO array at a spatial resolution of 
~ 11", Beuther et al. 2002b), we now attempted to derive 
radial intensity profiles from these higher spatial resolution 
data. Since the beam is elliptical, we extracted a strip in 
the east-west direction (smallest beam size of 0.5") and 
derived the profiles in the directions of positive and neg- 
ative R.A. offsets, respectively. The results and some fits 
to the data are shown in Figure 3. Obviously, the inten- 
sity profiles do vary significantly in the two directions with 
power-law indices between —1.0 and —3.6. Thus, the typ- 
ical approach of circular averaging is not useful here. Nev- 
ertheless, to better quantify the errors often introduced by 
circular averaging, we restored the image with a circular 
beam of ~ 0.75" and derived the intensity profile using cir- 
cular averaging. Not surprisingly, the derived power law 
index (not shown) is between the ones shown in Figure 3 
with an index of —2.1. One result to be kept in mind is 
that circular averaging has to be done cautiously, at least 
at this spatial resolution. 

Table 1 shows the submm continuum peak .Speak and in- 
tegrated intensity Si„ t . Comparing these values with the 
previously derived continuum flux densities (Bcuther et al. 
2004a), we find nearly the same peak intensity but about 
three times the integrated flux densities. The synthesized 
beam of the new dataset is a bit smaller than that of last 
years observations (l.l"x0.5" vs. 1.4"x0.9"). However, as 
shown in Fig. 2, data > 150 kA do not contribute strongly 
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to the measurable flux densities, and thus finding similar 
peak fluxes densities in both datasets is reassuring that 
the calibration of the SMA is reliable. The fact that we 
observe significantly higher integrated flux densities in the 
new dataset is due to the better uv-coverage we have this 
time. Due to this better uv-coverage we suffer less from 
missing short spacings and recover more of the extended 
halo emission. Nevertheless, measuring the integrated flux 
densities of the 850 /im single-dish SCUBA observations 
(Williams et al. 2004) in an aperture with a radius of 
~ 30", we get an integrated flux density of ~ 30 Jy. Thus, 
we only recover about 20% of the single-dish flux densi- 
ties with the SMA, the remaining 80% mostly reside in a 
larger-scale halo filtered out by the interferometer, and a 
smaller fraction of the integrated flux is lost because of the 
negative side-lobe emission. 

Following Hildcbrand (1983), assuming the submm con- 
tinuum emission to be due to optically thin dust emission, 
we calculate the gas mass and column density of the central 
core. The procedure and equations we use arc outlined in 
Bcuther et al. (2002b). Based on the spectral index anal- 
ysis from Beuther et al. (2004a) we assume a dust opacity 
index (3—1. Furthermore, we calculate the gas mass M 
on the one hand using a typical hot core temperature of 
100 K to better compare it with the results by Beuther 
et al. (2004a), and on the other hand using the tempera- 
ture of 350 K derived from the CH3CN observations (see 
§3.2.1). As discussed in (Bcuther et al. 2002b), the errors 
are dominated by systematics, e.g., the exact knowledge 
of the dust opacity index [3 or the temperature. While we 
take different temperatures into account the other param- 
eters are more difficult to acknowledge and we estimate 
the masses and column densities to be correct within a 
factor 5. 

Assuming 100 K we get a gas mass of 45 M Q close to the 
38 M© derived last year from the 1.3 mm continuum ob- 
servations (Beuther et al. 2004a), whereas using the 350 K 
results in a lower gas mas of 12 M©. The values will be 
compared below with the masses derived from the kine- 
matic HCOOCH3 observations (see §3.2.2). Regarding the 
gas column density N , both temperatures result in values 
of the order 10 cm~ 2 , close to the column densities previ- 
ously derived from the lower spatial resolution 1.3 mm data 
(Bcuther et al. 2004a). Assuming a relation between col- 
umn density and extinction of A„(mag) = N/(10 21 cm~ 2 ) 
we get a visual extinction of the order 1000. 

3.2. Spectral line data 

Figure 4 shows vector-averaged spectra of the com- 
plete spectral setup taken on a short baseline of 26 m, 
and Table 2 lists the detected lines. Line identification 
were conducted with the 350 GHz spectral line survey by 
Schilke et al. (1997a) and the molecular line databases 
JPL (Poynter & Pickett 1985) and CDMS (Muller et al. 
2001). Altogether 50 spectral lines from 18 species are 
observed (we consider different isotopomers as different 
species: CH 3 OH, 13 CH 3 OH, S0 2 , 34 S0 2 , SO, 34 SO, 
C2H5OH, CH3OCH3, NH 2 CHO, C 33 S, CN, OCS, H 2 CS, 
HCOOCH3, CH 3 CN, CTLf CN, CCH, HNCO). Interest- 
ingly, within this single spectral setup we detect molecular 
lines spanning upper state energy levels E uppcr from 17 to 
747 K (Table 2). Fourteen species are sufficiently strong to 



be imaged. Figure 5 shows the integrated images of these 
species, and only HCOOCH3 does peak exactly on the con- 
tinuum image. A more detailed analysis of the HCOOCH3 
data is given in §3.2.2. The lines we were originally mainly 
interested in, the CH3CN(19fe — 18k) K-ladder, appear to 
be optically thick, and do not properly trace the most cen- 
tral core (see §3.2.1). 

Comparing the other molecular lines with the previ- 
ously observed SiO outflow emission leaving the core in 
the northern direction (Beuther et al. 2004a), some lines 
show contributions from the outflow as well, most strongly 
H 2 CS. Most of the other species appear to peak close to 
the continuum emission - usually within < 1" correspond- 
ing to < 3600 AU - but none directly toward it. Interest- 
ingly, for most species this offset is toward the west and 
north-west, similar to the peculiar offset of the SiO outflow 
axis (Bcuther et al. 2004a). Since we do not detect any 
companion source in the submm continuum data (see also 
§4.1), it is unlikely that a secondary source is causing this 
offset. It appears more plausible that shock interactions 
of the outflow with the ambient molecular core enhance 
the molecular abundances and thus at least partly cause 
the observed offsets. It is difficult to quantify the relative 
contributions from excitation, chemistry, opacity, outflow 
and other kinematic effects to the morphological differ- 
ences between the various species, and likely an interplay 
of these effects has to be taken into account to explain all 
observational features. 

Comparing the spatial extent of the 50% contour lev- 
els from some of the molecules in Table 3 (we selected 
just the strong ones with reasonable S/N ratio), we find 
that HCOOCH3 is confined most strongly, approximately 
within the same region as the 860 [im continuum emis- 
sion, whereas all other molecules, even density tracers 
like CH 3 CN/CH 3 OH/C 33 S, are more extended. With 
HCOOCH3 peaking exactly toward the submm continuum 
peak and showing the smallest spatial extend, it appears 
to be one of the best dense gas tracers within this mini- 
line-survey. 

3.2.1. C'H 3 CN{19 k - 18 k ) 

All CH 3 CN(19fe — 18k) velocity channels peak offset from 
the mm continuum peak, and we cannot distinguish a clear 
velocity structure (Fig. 6). The picture is similar for all 
8 detected CH 3 CN(19 fc - 18 fe ) K-lincs (K=0,...,7). While 
physical reasons like the above mentioned outflow contri- 
butions can account for some of the offsets, the rather 
random-like distribution of offsets in this channel map is 
indicative of high optical depths throughout the whole K- 
ladder (see also the discussion about the spectral fits be- 
low). Thus, the CH 3 CN(19fe — 18&) lines appear not to be 
suited well as a disk tracer for this source. 

Nevertheless, we can use the CH 3 CN K-ladder to get 
a temperature estimate of the core. For this purpose, 
we used the XCLASS superset of the CLASS program 
(Schilke, priv. comm.). This software-package uses the line 
catalogs from JPL and CDMS, Poynter & Pickett 1985; 
Muller et al. 2001). We tried to fit the whole CH 3 CN 
spectrum simultaneously assuming optically thin emission 
in local thermodynamic equilibrium. Figure 7 presents 
the whole CH 3 CN(19fc — 18fc) spectrum toward the submm 
peak position and a model spectrum using a temperature 
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observations (Bcuther et al. 2004a). With the rotation 
axis, defined by the high-density tracer HCOOCH 3 , being 
approximately perpendicular to the outflow axis, these ob- 
servations in IRAS 18089-1732 show a disk signature sim- 
ilar to those known from many low-mass sources (e.g., D 
Tauri, Guillotcau & Dutrey 1998). This confirms at higher 
spatial resolution our previously proposed disk scenario for 
this High-Mass Protostellar Object. While the HCOOCH3 
observations at 216.966 GHz showed a P.A. of <~ 55° (from 
north), the axis of the new data is closer to 90°. Consider- 
ing the different energy levels of the various HCOOOCH3 
transitions (£ U p per (HCOOCH 3 (202 i 9 -19 19)) ~ 111K & 
£ U pper(HCOOCH 3 (289,20 - 27 9 ,i 9 ) j ~ 294 K), this implies 
that the lower-energy-level transitions are more strongly 
affected by the infalling surrounding envelope. In contrast, 
the higher-energy-lines appear to be better probes of the 
actual underlying disk properties. 

The deconvolved source size of the integrated 
HCOOCH3 map in Fig. 5 is - 2050 AU (using the 50% 
contour as the observed image size), similar to the separa- 
tion of the blue- and red-shifted HCOOCH3 peaks in Fig. 9 
of ~ 1800 AU. With this source size, we cannot clearly 
distinguish whether we observe the outer regions of a cir- 
cumstellar disk or rather a larger-scale rotating molecular 
envelope/torus which could feed the inner accretion disk. 
Assuming equilibrium between the centrifugal and gravi- 
tational forces at the outer radius of the disk the enclosed 
dynamical mass can be estimated via: 
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of 350 K and a CH3CN column density of 1.4 x 10 15 cm~ 2 . 
This model fits the data good for the K=0,l,2,4,5,7 lines 
but bad for K=3,6. The ortho-CH 3 CN (A-species) lines 
K=3,6 have double the statistical weight than the other 
observed lines, whereas the observed line strengths of the 
K=3,4,5,6 lines (including para and ortho species) are sim- 
ilar. The observed line ratio close to unity between ortho- 
and para-CH 3 CN lines is additional evidence for the large 
optical depths of the lines. In spite of these difficulties, 
the temperature of 350 K is consistent with our previous 
temperature estimate based on HCOOCH3 observations 
(Bcuther et al. 2004b) and thus gives a reasonable tem- 
perature for the central forming hot molecular core. 

3.2.2. Disk signature in HCOOCH 3 

In contrast to the CH3CN emission, HCOOCH3 traces 
the central dust emission well and thus appears to be 
likely optically thin. Unfortunately, the spectral line 
at ~ 348.91GHz is weaker than most of the CH 3 CN 
lines, and it is a blend of two HCOOCH3 lines and the 
CH 3 CN(19 9 - 18 9 ) line (see Table 2). However, the 
energy levels of the HCOOCH3 lines are at ~ 294 K 
whereas for CH 3 CN(19 9 - I89) it is at - 745 K (Ta- 
ble 2). Therefore, the spectrum is dominated by the 
HCOOCH3 lines. Furthermore, using the rest-frequency of 
the HCOOCH 3 (289,2o-27 9 ,i9)-E line at 348.9095 GHz, the 
line center is at 32km/s consistent with the line centers of 
other dense gas tracers, e.g., CH 3 CN. The frequency of the 
corresponding HCOOCH 3 (28g,2o — 27g,ig)-A line should be 
blue-shifted by ~ 4.7km/s but we do not find a signif- 
icant secondary peak toward this velocity. Likely, there 
is underlying emission from this A-type line, although at 
weaker intensities, but the spectrum is dominated by the 
HCOOCH3 line of the corresponding E-species. Thus, it 
appears feasible to use this HCOOCH3 line for a closer 
analysis of the kinematics of the central core. 

Figure 8 shows a channel map of this 
HCOOCH3 (289,20 - 27 9 ,i 9 )-E line. Consistent with our 
previous investigation (Beuthcr et al. 2004a), we find a 
velocity shift in the east-west direction over the core 3 
. This can be depicted in the blue channels at 26 and 
28km/s, which show compact emission close to and a 
bit east of the mm continuum peak, and the red chan- 
nel at 34km/s, which shows compact emission west of 
the submm continuum peak. The intermediate chan- 
nels at 30 and 32km/s show slightly more extended 
emission features. The fact that we see blue emission 
over a little bit wider velocity range than red emission 
may be due to the weaker HCOOCH3 (289,20 — 27g,i 9 )- 
A line, which is blue-shifted with regard to the main 
HCOOCH 3 (28 9 ,2o - 27 9 ,i 9 )-E line. The red emission of 
the HCOOCH 3 (28 9 ,2o - 27 9 ,i 9 )-A line would then be un- 
detectable because of the stronger ambient gas emission 
at 30 and 32km/s from the HCOOCH 3 (289, 2 o - 27 9 ,i 9 )-E 
line. Figure 9 presents an overlay of the HCOOCH3 data 
binned into blue/central/red components. Additionally, 
the outflow axis is sketched as traced by the SiO(5-4) 

3 We note that one cannot directly compare the velocity channels of the previous HCOOCH3 observation near 216.966 GHz with our new data, 
because in last year's analysis we set the rest-frequency to a central component of the line blend whereas for the new data it is set to the 
strongest component HCOOCH3(289,2o — 27g,ig)-E. This introduces slightly different velocity ranges for both observations, but it does not 
change the relative velocity shift over the source. 

4 The HWZI of the whole observed line is ~ 5km/s. Since there is some blending between the HCOOCH3(28g,20 — 27g,ig) main E-type and 
weaker A-type line, we adopt a slightly smaller value of HWZI=~ 4km/s for the HCOOCH3 (289,20 — 27g,ig)-E transition. 
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where Mdi sk is the disk mass, r the disk radius, in our 
case half the separation of the blue and red HCOOCH3 
emission peaks at most extreme blue and red velocities, 
and v the Half Width Zero Intensity (HWZI) of the spec- 
tral line. Equations 1 & 2 have to be divided by sin 2 (i) 
where i is the unknown inclination angle between the 
disk plane and the line of sight. Adopting the values 
v ~ 4 km/s and r ~ 900 AU, the resulting disk mass 
^disk is ~ 16/(sin 2 (i)) M© 4 . This value is close to the 
22/ (sin 2 (i)) M© value derived from our previous lower res- 
olution 1.3 mm HCOOCH3 observations (Beuther et al. 
2002c). It also compares well to the 12 M core mass de- 
rived assuming 350 K (in contrast to the estimated 45 M© 
assuming 100 K, see §3.1 and Table 1). As already noticed 
by several studies (e.g., Cesaroni et al. 1997; Zhang et al. 
1998; Sandell et al. 2003; Beltran et al. 2004; Beuther et al. 
2004a) , in contrast to low-mass disks where the disk mass 
is usually negligible compared with the core mass, a con- 
siderable amount of the gas mass appears to be associated 
with large accretion disks and/or rotating envelopes in the 
central regions of high-mass star formation. 

3.3. Varying spatial distributions of CH3OH 
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The previously observed CH30H(5i,4 — 42,2) line at 
216.946 GHz peaked toward the mm continuum and exhib- 
ited a velocity gradient in the east-west direction (Fig. 3 
in Bcuther et al. 2004b). Figure 10 presents channel 
maps for the CH 3 OH(2 2 , - 3i, 3 ), CH 3 OH(7 5 , 3 - 64,2), 
and CH 3 OH(14 M3 - 14 ,i 4 ) lines at 340.141, 338.722, and 
349.107 GHz, respectively. The spatial distribution of the 
three lines is not exactly the same: the CH 3 OH(2 2 ,o — 3i ;3 ) 
and CH 3 OH(14i,i 3 — 14o.i4) both show red emission to 
the west of the continuum peak but only weak or barely 
detectable blue emission to the east (Fig. 10 top and bot- 
tom). Contrary to this we cannot discern such a signature 
in the CH 3 OH(75 j3 — 64,2) line, the emission peaks north 
and south of the submm continuum peak and seems to be 
affected by the outflow in that direction (Fig. 10 middle). 
While none of the lines shows a clear east- west velocity gra- 
dient indicative of the proposed disk as suggested by the 
previous CH 3 OH(5i 4 — 4 2 2) observations, the red emis- 
sion of the CH 3 OH(2 2 , - 3i', 3 ) and CH 3 OH(14 M3 - 14 ,i 4 ) 
lines are at least consistent with it. 

Since CH 3 OH is a slightly asymmetric rotor with a very 
rich spectrum in the mm and submm bands (see Table 2 or 
Leurini et al. 2004 for a more detailed CH 3 OH discussion), 
it is interesting to investigate in what kind of environment 
one finds the molecule. As shown in Table 2, within the 
given bandpass we observe CH 3 OH lines spanning upper 
level energy ranges between 44 and 695 K. Since the lower 
energy CH 3 OH(2 2 , -3i, 3 ) (£ upP cr = 44 K) and the higher 
energy CH 3 OH(14i j i 3 - 14 ,i 4 ) (E uppcr = 259 K) lines are 
more confined to the inner center than the intermediate 
energy CH 3 OH(75 i3 — 64^) line (-E U ppcr = 87 K), one could 
speculate that the temperatures within the outflow might 
be in the latter regime and thus influence this line most 
strongly. However, not all lines of the CH 3 OH(7 — 6) series 
resemble the outflow-influenced morphology but others 
also trace the inner core again (e.g., CH 3 OH(74 i3 — 64, 3 )). 
In the literature one finds cases where CH 3 OH traces the 
protostellar cores (e.g., Wyrowski et al. 1999; van der Tak 
et al. 2000) as well as molecular outflows (e.g., Bachillcr 
et al. 1998; Beuther et al. 2002a; J0rgensen et al. 2004), 
and currently it is not well understood which CH 3 OH line 
is suitable to trace one or the other effect. Nevertheless, 
because CH 3 OH is a molecule with often more than one 
line in a given spectral setup, it is a potentially useful tool 
to investigate the different processes. Since the CH 3 OH 
collisional rates were calculated recently (Pottage et al. 
2004), the potential of CH 3 OH as an interstellar tracer is 
likely going to be exploited more extensively than before 
(e.g., Leurini et al. 2004). 

4. DISCUSSION 

4.1. Continuum emission 

It is interesting to see that even at this high spatial 
resolution IRAS 18089-1732 still does not split up in any 
sub-fragments. The 3<r continuum sensitivity of 105 mJy 
corresponds to a mass sensitivity of 0.3, 0.7, or 3.0 M Q , 
depending on the assumed temperature which we set to 
350 K (the derived CH 3 CN hot core temperature), 100 K 
(a more typical hot core temperature), and 30 K (which 
should be the temperature of the surrounding envelope), 

5 It is noteworthy that the mm continuum peak in IRAS 20126+4104 
1999). This may be responsible for the large-scale precession of the mt 



respectively. Since potential lower-mass companions of an 
evolving massive cluster should be distributed a little bit 
further out from the central core (e.g., Lada & Lada 2003), 
lower temperature values of 100 and/or 30 K appear to be 
a more realistic for potential companions. Nevertheless, 
we find no companion down to a 3cr mass limit of 0.7 and 
3M Q at 100 and 30 K, respectively. The non-detection 
of any companion source supports the previous reasoning 
that the SiO outflow is indeed emanating from the main 
massive submm continuum source and not caused by a 
lower-mass cluster member (Beuther et al. 2004a). The 
small offset of the SiO outflow axis from the submm peak 
position (~ 1", Beuther et al. 2004a) is still a bit puz- 
zling, but it can likely be attributed to the underlying 
flow geometry and the resulting shock-induced local SiO 
enhancements (e.g., Schilke et al. 1997b). 

The single submm continuum peak is also distinctively 
different to the recent study by Beuther & Schilke (2004) 
who found the HMPO IRAS 19410+2336 to split up into a 
cluster of sources with a protocluster mass function resem- 
bling the stellar Initial Mass Function. Contrary to this, 
there exist other prominent examples which do not split up 
into clusters, e.g., IRAS 20126+4104, but remain one mm 
continuum source down to the highest spatial resolution 
(Cesaroni et al. 1999) 5 , similar to the case of IRAS 18089- 
1732. 

What is causing these observational differences? On 
the one hand, it is possible that the study by Beuther 
& Schilke (2004) had very favorable observational circum- 
stances, because they observed IRAS 19410+2336 in three 
configurations with the Plateau de Bure Interferometer 
and thus had an exceptionally good uv-coverage sampling 
various spatial scales. Our SMA data have been observed 
with baselines not shorter than 20 m and thus does not 
trace the corresponding larger spatial scales, especially at 
this high frequencies. Excluding the compact Plateau de 
Bure configuration, Beuther & Schilke (2004) would not 
have detected that many sub-sources in IRAS 19410+2336 
as well. Therefore, observational differences might ac- 
count to some degree for the discussed differences. On 
the other hand, could it be possible that some massive 
star-forming regions evolve in a fragmenting cluster-like 
way like IRAS 19410+2336, whereas other star-forming re- 
gions exhibit the core- halo structure like IRAS 18089-1732 
at the earliest evolutionary stages and develop the cluster- 
like structure at later evolutionary stages? This scenario 
would imply that at least for such kind of sources the mas- 
sive cluster members would form first and the low- and 
intermediate-mass members later. One way to solve this 
question is to choose a number of similar sources (similar 
in the sense of evolutionary stage, luminosity and distance; 
IRAS 20126+4104 and IRAS 19410+2336 for example ful- 
fill these criteria) and observe all these sources with ex- 
actly the same intcrfcrometric configurations to trace the 
same spatial scales. If this experiment would still result 
in the same observational differences as encountered here, 
this would be a hint to the existence of a variety of evolu- 
tionary paths for massive star-forming regions. 

Another result from the continuum study is that it is dif- 
ficult to derive a consistent intensity distribution and thus 

contains a ionized double jet indicating a binary system (Hofncr et al. 
lecular outflow (Shepherd et al. 2000; Cesaroni et al. 2005). 
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density distribution on this small spatial scales because 
apparently IRAS 18089-1732 is not symmetric at all. It 
is likely that the central core still harbors multiple unre- 
solved sources which could cause these asymmetries. In 
any regard, the usual approach of circular averaging and 
then deriving density distributions which are interpreted 
within one or the other theoretical model has to be ap- 
proached cautiously and needs to be re-checked for each 
studied source separately. 

4.2. Disks in massive star formation 

We are able to confirm the previously suggested 
disk/outflow scenario for the HMPO IRAS 18089-1732 at 
higher spatial resolution. The HCOOCH3 velocity shift 
is approximately perpendicular to the SiO outflow axis. 
However, even at this spatial resolution it is hard to judge 
whether we really trace a Keplerian accretion disk (Keple- 
rian in the sense of a Keplerian velocity profile dominated 
by the central protostar) or rather a rotating larger-scale 
envelope which rotates perpendicular to the outflow but 
which might not be Keplerian at all. The fact that a large 
fraction of the central gas core appears to be associated 
with the rotation indicates that it might not be in Keple- 
rian motion. This implies a follow-up question whether 
massive accretion disks need to be Keplerian like their 
low-mass counterparts, or whether massive disks may be 
non-Keplerian, potentially self-gravitating structures (self- 
gravitating in the sense that the gravitational potential of 
the disk-star system is dominated by the disk). While the 
observations toward prototypical high-mass disk source 
IRAS 20126+4104 favor a Keplerian disk, although with 
a central mass of ~ 7M it is not a very massive source 
(e.g., Cesaroni et al. 2005), other sources indicate disk 
masses well in excess of the central source (e.g., Sandcll 
et al. 2003; Chini et al. 2004). For IRAS 18089-1732, we 
find a rotating structure which is consistent with an ac- 
cretion disk/rotating envelope, but unfortunately the data 
still lack the spatial resolution to study the physical struc- 
ture of the disk. Follow-up observations at higher angular 
resolution are one tool to tackle these questions. 

In addition, the question posed in §1 to find the most 
suitable disk-tracing molecules has not been sufficiently 
answered. The previously successful molecule CH3CN has 
proved in this source as not suitable. It is not clear whether 
simply the chosen K-laddcr line scries is a bad one because 
the opacity r scales with the frequency v like: 

<^> r = constant x v — lj N u (4) 

with A the Einstein coefficient, Av the line width, and 
N u the column density of the upper state. Using the tem- 
perature estimate of 350 K and assuming the same line 
width for all CH3CN line series, we can approximate the 
change in optical depth for the CH3CN line series in the 
349 GHz and 220 GHz bands by: 

T Cg 3 CjV(19-18) ^ 1 q x x ^u,CH 3 CN(19-!S) 

1~CH 3 CN(12-11) 220 N UtCH3CN (12-11) 

In case the CH 3 CN column densities in both transitions 
are similar, which is given to first order in thermal equilib- 



rium at the estimated temperature of 350 K, the line opac- 
ity in the 860 um band would be approximately 2.5 times 
larger than in the 1 mm band. Therefore, higher optical 
depth at 860 /im could contribute to the observation that 
we do not find a disk signature in the CH3CN(19fc — 18k) 
transitions. However, it is possible that already the 1mm 
lines are optically thick which has to be checked via future 
observations of the latter line series. 

Contrary to this, the HCOOH3 molecule traces the disk- 
like structure in the 1mm and the 860 um band success- 
fully and apparently does not suffer from too high opacity. 
However, in spite of its ubiquity in many spectral bands 
(Sutton et al. 1985; Schilke et al. 1997a, 2001), spectral 
HCOOCH3 lines are often rather weak and and blended 
with various components, and therefore observationally 
more difficult to tackle. To start statistically significant 
studies of massive accretion disks, we still have to find 
a suitable molecular transition which unambiguously can 
trace disks and does not require a too large amount of 
observing time. 

In spite of these problems, it is promising that more and 
more massive disk candidates are being studied. With con- 
certed efforts from various groups employing all available 
interferometers, it is likely that we will get a much better 
understanding of massive disks over the next decade. 



5. CONCLUSIONS AND SUMMARY 

Fifty spectral lines from eighteen different species span- 
ning an energy range from 17 to 747 K were detected 
within the 4 GHz bandpass at 860 /im. While the pre- 
viously assumed best disk tracer CH3CN in this band is 
optically thick and thus not suitable for follow-up disk 
studies in this source, HCOOCH3 appears optically thin 
and shows a velocity-disk signature similar to well-known 
low-mass disk sources. The estimated mass of this rotating 
structure, assuming equilibrium between the gravitational 
and centrifugal forces, is 16/(sin 2 (i)) M©, of the same or- 
der as the core mass derived from the submm dust con- 
tinuum emission. This implies that a significant fraction 
of the central core mass is associated with the rotation, 
and the assumed disk-like structure may not be Keple- 
rian. Within this mini-line-survey, HCOOCH3 appears as 
one of the best disk-tracing molecules. 

Although CH3CN(19fe — 18k) is optically thick, we can 
use the K-ladder to estimate a temperature in the region. 
We find a temperature of approximately 350 K, consistent 
with previous HCOOCH3 estimates (Beuther et al. 2004b), 
and a hot molecular core has already formed. 

Even at 1" spatial resolution, we do not detect any com- 
panion source in the field of view down to mass limits 
between 0.2 and 3M Q (depending on the assumed tem- 
perature). Since this finding is distinctively different to 
the recent study toward the HMPO IRAS 19410+2336 it 
might indicate different evolutionary paths for the forma- 
tion of massive clusters. However, observational systemat- 
ics can also account for some of the differences, and more 
detailed investigations of various sources with the same ob- 
servational settings are necessary to draw firm conclusions. 
Furthermore, the central dust core is highly non-circular, 
indicating potential difficulties for circular averaged stud- 
ies of intensity profiles at this high spatial resolution. 



SMA sub-mm observations of IRAS 18089-1732 



7 



We like to thank Todd Hunter for the fruitful discus- 
sions about this source. Thanks a lot also to Peter Schilkc 
for providing the XCLASS software. Furthermore, we ap- 
preciate the constructive and quick report from an anony- 



mous referee. H.B. acknowledges financial support by the 
Emmy-Noether-Program of the Deutsche Forschungsge- 
meinschaft (DFG, grant BE2578/1). 



REFERENCES 



Bachiller, R., Codclla, C, Colomer, F., Liechti, S., & Walmsley, C. M. 

1998, A&A, 335, 266 
Beltran, M. T., Cesaroni, R., Neri, R., et al. 2004, ApJ, 601, L187 
Bcuther, H., Hunter, T. R., Zhang, Q., et al. 2004a, ApJ, 616, L23 
Bcuther, H. & Schilke, P. 2004, Science, 303, 1167 
Beuther, H., Schilke, P., Gueth, F., et al. 2002a, A&A, 387, 931 
Bcuther, H., Schilke, P., Menten, K. M., ct al. 2002b, ApJ, 566, 945 
Bcuther, H., Schilkc, P., Sridharan, T. K., et al. 2002c, A&A, 383, 

892 

Beuther, H., Walsh, A., Schilke, P., et al. 2002d, A&A, 390, 289 
Bcuther, H., Zhang, Q., Hunter, T. R., et al. 2004b, ApJ, 616, L19 
Cesaroni, R., Felli, M., Jenness, T., et al. 1999, A&A, 345, 949 
Cesaroni, R., Felli, M., Testi, L., Walmsley, C. M., & Olmi, L. 1997, 
A&A, 325, 725 

Cesaroni, R., Neri, R., Olmi, L., Walmsley, C, & Hofner, P. 2005, 
A&A accepted 

Chini, R., Hoffmeister, V., Kimeswenger, S., et al. 2004, Nature, 429, 
155 

Davis, C. J., Varricatt, W. P., Todd, S. P., & Ramsay Howat, S. K. 

2004, A&A, 425, 981 
Garay, G., Faiindez, S., Mardoncs, D., ct al. 2004, ApJ, 610, 313 
Gibb, A. G., Hoare, M. G., Little, L. T., & Wright, M. C. H. 2003, 

MNRAS, 339, 1011 
Guilloteau, S. & Dutrey, A. 1998, A&A, 339, 467 
Hildebrand, R. H. 1983, QJRAS, 24, 267 

Hofner, P., Cesaroni, R., Rodriguez, L. F., & Marti, J. 1999, A&A, 
345, L43 

j0rgensen, J. K., Hogerheijde, M. R., Blake, G. A., et al. 2004, A&A, 
415, 1021 

Lada, C. J. & Lada, E. A. 2003, ARA&A, 41, 57 

Leurini, S., Schilke, P., Menten, K. M., et al. 2004, A&A, 422, 573 



Miiller, H. S. P., Thorwirth, S., Roth, D. A., & Winnewisser, G. 

2001, A&A, 370, L49 
Pestalozzi, M. R., Elitzur, M., Conway, J. E., & Booth, R. S. 2004, 

ApJ, 603, L113 

Pottage, J. T., Flower, D. R., & Davis, S. L. 2004, MNRAS, 352, 39 
Poynter, R. L. & Pickett, H. M. 1985, Appl. Opt., 24, 2235 
Sandell, G., Wright, M., & Forster, J. R. 2003, ApJ, 590, L45 
Schilke, P., Benford, D. J., Hunter, T. R., Lis, D. C, & Phillips, 

T. G. 2001, ApJS, 132, 281 
Schilke, P., Groesbeck, T. D., Blake, G. A., & Phillips, T. G. 1997a, 

ApJS, 108, 301 

Schilkc, P., Walmsley, C. M., Pineau des Forets, G., & Flower, D. R. 

1997b, A&A, 321, 293 
Scoville, N. Z., Carlstrom, J. E., Chandler, C. J., et al. 1993, PASP, 

105, 1482 

Shepherd, D. S., Yu, K. C, Bally, J., & Testi, L. 2000, ApJ, 535, 833 
Sridharan, T. K., Beuther, H., Schilke, P., Menten, K. M., & 

Wyrowski, F. 2002, ApJ, 566, 931 
Su, Y.-N., Liu, S.-Y., Lim, J., et al. 2004, ApJ, 616, L39 
Sutton, E. C, Blake, G. A., Masson, C. R., & Phillips, T. G. 1985, 

ApJS, 58, 341 

Torrelles, J. M., Gomez, J. F., Rodriguez, L. F., ct al. 1996, ApJ, 
457, L107+ 

van der Tak, F. F. S., van Dishoeck, E. F., & Casclli, P. 2000, A&A, 
361 327 

Williams, S. J., Fuller, G. A., & Sridharan, T. K. 2004, A&A, 417, 
115 

Wyrowski, F., Schilkc, P., Walmsley, C. M., & Menten, K. M. 1999, 
ApJ, 514, L43 

Zhang, Q., Hunter, T. R., & Sridharan, T. K. 1998, ApJ, 505, L151 
Zhang, Q., Hunter, T. R., Sridharan, T. K., & Ho, P. T. P. 2002, 
ApJ, 566, 982 



2 - 3600AU 




♦ 

_4 I i I i I i i i I i i i I i L_ 

4 2 -2 

Aa ["] 

Fig. 1. — 860 /im continuum image of IRAS 18089-1732. The contours range from 175 to 1400 mjy in 175 mjy steps. The beam (1.1" X 
0.5" P.A. 1°) is shown at the bottom-left. The deconvolved source size is ~ 3000 AU (using the 50% contour as the observed image size). The 
axis offsets are in arcseconds from the phase center (§2). 
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Fig. 2. — Amplitude versus uv-distance plots of the continuum data presented in Figure 1. The left panel shows the whole dataset in 5kA 
bins, the right panel presents only the long baseline data in larger 40 kA bins. The dotted lines show the expected amplitude for zero signal 
with the given statistical errors. 




Fig. 3. — Continuum flux versus peak-offset plot. The data arc taken along a strip through the continuum peak in east-west direction - 
the direction with the smallest beam-size of 0.5". The filled triangles mark positive offsets and the open squares negative offsets. The dashed 
lines show power-law fits to subsets of the data, outlining that there is significant variation even along this singular strip. 



Table 1 

sltbmillimeter continuum data 



*5pcak 


1.400 mJy 




6440 mJy 


M at 100 K 


45 M Q 


N at 100 K 


3.5 x 10 23 cm~ 2 


M at 350 K 


12 M 


N at 350 K 


9.5 x 10 23 cm- 2 



SMA sub-mm observations of IRAS 18089-1732 
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Fig. 4. — Vector-averaged spectra on a short baseline of 26 m. Most molecular lines from Table 2 are marked. 
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Fig. 5. — Integrated images of all observed molecular species. The contour levels are chosen for each image from ±10 to ±90% of its peak 
emission. The star marks the submm continuum peak position (Fig. 1). 
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Fig. 6. — Channel map of CH3CN(192 — I82). The grey-scale shows the 860 fim continuum image and the contours present the CH3CN 
data. Positive contour levels (full lines) are 0.67 to 4.69 Jy (0.67 Jy steps), and negative contours range from —0.67 to —2.01 Jy also in —0.67 Jy 
steps. 
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Fig. 7. — The full line shows the CH3CN(19fc — 18&) spectrum taken toward the submm continuum peak after imaging the lines. The 
dotted line presents the model spectrum at 350 K with a column density of 1.4 X 10 15 cm -2 . The K-levels are indicated above each line. 




Aa ["] 

Fig. 8. — Channel map of HCOOCH3. The positive contour levels are in solid lines from 0.35 to 1.75 Jy (0.35 Jy steps), the negative 
contours are in dashed lines from -0.35 to -1.05 Jy (0.35 Jy steps) The central velocities are marked at the top-right of each panel, and the 
star pinpoints the submm continuum position (Fig. 1). 
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Fig. 9. — Disk signature in HCOOCH3. Dashed and dotted lines show the blue and red disk emission, respectively. The full lines present the 
emission at intermediate velocities. Again, the contour levels are chosen from 10 to 90% of the peak emission of each velocity bin separately. 
The bold lines marks the outflow direction observed in SiO(5-4) (Beuthcr et al. 2004a). 
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Fig. 10.— Channel maps of the CH 3 OH(2 2 ,o - 3i, 3 ) (top), CH 3 OH(7 5 ,3 - 64,2) (middle), and CH 3 OH(14i,i 3 - 14 ,i4) (bottom) lines, 
respectively. The grey-scale shows the 860 fim continuum emission in all panels. The CH 3 OH(22,o — 3i j3 ) map is contoured positive from 0.46 
to 3.68 Jy (step 0.46 Jy), and negative from -0.46 to -1.84 Jy (step -0.46 Jy). The CH 3 OH(75. 3 — 64,2) date are contoured positive from 0.6 
to 4.8 Jy in 0.6 Jy steps, and negative at —0.6 & — 1.2 Jy. The CH 3 OH(14i i i 3 — 14o,i4) arc contoured positive from 0.76 to 5.32 Jy in 0.76 Jy 
steps and negative at —0.76 and — 1.52 Jy. The synthesized beams are shown at the bottom-left of each panel (full ellipse for continuum and 
lined ellipse for CH 3 OH). 
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Table 2 



Observed lines 



V 


line 


Ob 

peak 

[K] 


E 


GHz 




[K] 


LSB 


338.5041 


CH 3 OH(7 2 , 6 - 6i, 5 )E 


5.0 


153 


338.5126 


CH 3 OH(7 6 , 2 - 6 5 .i)A 


4.1 


145 


338.5303 


CH 3 OH(7 6;2 - 6 5 ,i)E 


3.7 


160 


338.5408 


CH 3 OH(7 5i2 - 6 5 . 2 )A+ a 


4.9 


115 


338.5432 


CH 3 OH(7 5 . 2 - 6 5 . 2 )A- a 




115 


338.5599 


CH 3 OH(7 2 , 5 - 6 2 , 5 )E 


5.8 


128 


338.5832 


CH 3 OH(7 5 , 2 - 6 5 , 2 )E 


5.3 


113 


338.6150 


CH 3 OH(7 4 , 3 - 6 4 , 3 )E 


5.0 


86 


338.6399 


CH 3 OH(7 5i3 - 6 4 , 2 )A+ 


4.9 


103 


338.7216 


CH 3 OH(7 5 , 3 - 6 4 , 2 )E a 


4.8 


87 


338.7229 


CH 3 OH(7 3! 5 - 6 2 , 4 )E a 




91 


338.7604 


id CH 3 OH(137 !7 - 12 7 , 6 )A+ 


2.2 


206 


338.7858 


34 O 1 A A 1 o \ 

M SO 2 (14 4) i - 12 3) i) 




134 


338.8862 


C 2 H 5 OH(15 7 , 8 - 15 6 , 9 ) a 




162 


338.8873 


C 2 H 5 OH(15 7 , 9 - 15 6 .io) a 




162 


339.3126 


C 2 H 5 OH(12 7 , 5 - 12 6;6 ) 




126 


339.3415 


SO(3 3 - 2 3 ) 


1.6 


25 


339.4917 


CH 3 OCH 3 (19i is - 18 2 i 7 )AA 


3.4 


176 


339.8576 


34 SO(8 9 - 7 8 ) 


2.6 


77 


339.9041 


NH 2 CHO(16 6ill - 15 6 ,io) 




246 


340.0527 


C 33 S(7-6) 


2.6 


65 


340.1412 


CH 3 OH(2 2 . - 3i, 3 )A+ 


4.4 


44 


340.1893 


C 2 H 5 OH(6 5 , 2 - 5 4 ,i) 


1.6 


49 


340.2486 


CN(3, 5 , 2 _ 2 ; 2' 2 ) 




32 


340.3164 


S02(282,26 — 281,27-) 


1.6 


392 


340.3937 


CH 3 OH(16 6 ,io - 17 5 .i 3 )A+ 


3.4 


508 


340.4492 


OCS(28 - 27) 2.3 


237 





USB 



348.5319 


H 2 CS(10i, 9 - 9i, 8 ) 


4.3 


105 


348.9095 


HCOOCH 3 (28 9 ,2o - 27 9 ,i 9 )E a 


1.7 


294 


348.9114 


CH 3 CN(19 9 - 18 9 ) a 




745 


348.9150 


HCOOCH 3 (28 9 .2o - 27 9 .i 9 )A a 




294 


349.0485 


HCOOCH 3 (28 9 .i 9 - 27 9 ,i 8 )E 




294 


349.1070 


CH 3 OH(14i,i 3 - 14 ,i4)A 


6.2 


259 


349.1253 


CH 3 CN(19 7 - 18 7 ) 


2.1 


517 


349.1732 


CH 3 13 CN(19 4 - 18 4 ) 




281 


349.2121 


CH 3 CN(19 6 - 18 6 ) 


3.8 


424 


349.2858 


CH 3 CN(19 5 - 18 5 ) 
CCH(4, |, | — 3, 2; 5) 
CH 3 CN(19 4 - I84) 


3.2 


346 


349.3379 




41 


349.3461 


3.4 


281 


349.3930 


CH 3 CN(19 3 - 18 3 ) 


3.4 


231 


349.4266 


CH 3 CN(19 2 - 18 2 ) 


4.8 


196 


349.4467 


CH 3 CN(19i - 18i) a 




175 


349.4534 


CH 3 CN(19 - 18 ) a 


6.1 


167 


350.1031 


13 CH 3 OH(l 1 , -0 ,o)A+ 




17 


350.2877 


CH 3 OH(15 3 ,i3 - 16 4 ,i 3 ) 


2.0 


695 


350.3333 


HNCO(16 M6 - 15i,i B ) 


3.0 


185 


350.4235 


CH 3 CN(w 8 = 1)(19_ 2 - 18_ 2 ) 


2.1 


747 


350.4423 


HCOOCH 3 (28 8 ,2o - 27 8 . 20 )E a 


2.1 


283 


350.4449 


CH 3 CN( V8 = 1)(19 - 18 ) a 




693 



a line-blcnd 

b The peak intensities Speak were measured toward the submm contin- 
uum peak position after imaging the whole data-cube with 2km/s spectral 
resolution. Lines which were too weak to image are marked with a 
and for line-blends we just list the intensity for one component. 
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Table 3 

Approximate area within 50% contours™ 



Area 
[arcsec 2 ] 



860 /im continuum 1.0 

HCOOCH 3 (289. 20 - 27 9 .i 9 )E 1.2 

CH 3 CN(19 2 - 18 2 ) 1.8 

CH 3 CN(w 8 = l)(19- 2 - 18_ 2 ) 1.6 

CH 3 OH(14 lil3 - 14o,i4)A 1.7 

13 CH 3 OH(13 7 7 - 12 7 6 )A+ 1.3 

CH 3 OCH 3 (19ii 8 -18 2 i 7 )AA 1.9 

H 2 CS(10i, 9 - 9i, 8 ) 2.7 

HNCO(16i.i 6 -15i.i 5 ) 1.6 

OCS(28 - 27) 3.8 

C 33 S(7-6) 2.3 



a The area is calculated assuming an elliptical 
source geometry and measuring the major and mi- 
nor axes of the 50% contour levels. 



